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A Method to Construct Asymptotic Solutions 
Invariant under the Renormalization Group 

Masatomo IWASA*) and Kazuhiro NOZAKI**) 
Departmant of Physics, Nagoya University, Nagoya 464-8602, Japan 

A renormalization group method with the Lie symmetry is presented for the singular 
perturbation problems. Asymptotic solutions are obtained as group- invariant solutions under 
approximate Lie groups admitted by perturbed diferential equations. 

§1. Introduction 

There have been many studies concerning application of the renormalization 
group theory of particle physics as a singular perturbation method to treat differen- 
tial equations since the work of the Illinois group. ^ In this method, integral constants 
appearing in the lowest-order perturbed solution are renormalized in order to remove 
secular or divergent terms appearing in the naive perturbation solution and give a 
well-behaved asymptotic solution. We call this renormalization group method in 
singular perturbation theory the "conventional RG method" in this paper. A geo- 
metrical aspect of the conventional RG method has been studied by means of the 
envelopes of solutions. 3 ) The conventional RG method and its variants depend more 
or less on the naive perturbation analysis. In connection to this, it is interesting 
that the geometrical symmetry of a differential equation with respect to continuous 
transformations, i.e. the Lie symmetry, is known to be useful for deriving new so- 
lutions from known solutions. 4 ) There have been some pioneering works aimed of 
constructing the renormalization group method in terms of the Lie symmetry. 5 ) 6 ) 
However, to this time, there is no case in which have been derived asymptotic solu- 
tions in singular perturbation problems by means of the Lie symmetry. The purpose 
of this paper is to derive such asymptotic solutions in the framework of the Lie group 
and symmetry. 

§2. Method of the renormalization group with symmetry 

For simplicity, let us consider the second-order ordinary differential equation 

F Q (u,u',u") + eF 1 (u,u') = 0, (2-1) 

where we have v! = ^ and u" = Fq and F\ are some functions of their argu- 
ments, and e is a small parameter. This system is assumed to be nearly solvable in 
the sense that the general solution is known for the unperturbed case, e = 0. Our 
purpose is to construct an approximate solution of the system H'iTjl for small e by 
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incorporating an approximate Lie symmetry. First, the equation (|2TJI is rewritten 
as a system of first-order ordinary differential equations in the form 

U> = U ' (2-2) 
F o (u,v,v') + eF 1 {u,v) = 0, v ' 

or as a first-order ordinary differential equation with respect to a complex valuable 
z := u + iv as 

G (z,z,z!,&)+eG 1 (z,z) = 0. (2-3) 

Let (|2-2[) admits a Lie group transformation whose infinitesimal generator takes the 
form 

X = 8 £ + C(t, u, v)8 t + rf{t, u, v)8 u + V v (t, u, v)d v , (2-4) 
Then, under the transformation (u, v) t—* (z,z), (|2-4[) can be rewritten as 

X = d e + £(t, z, z)d t + r] z (t, z, z)d z + rf{t, z, z)d- z , (2-5) 

where 

£ = £ = rf = rj u + ir] v , rf = rj u - irf = rf . (2-6) 

Next, let 

z, z, z', z) = 8 £ + £(t, z, z)d t + r] z (t, z, z)d z + r] z (t, z, z)d- z 

+r] z W (t, Z) z, z', z')d z , + n z ^ (t, z, z, z', z')d- z , (2-7) 

be the prolongation of X. Then, the determining equation for (j2-3j) . which determines 
each component of the vector field (|2-5(l . is given by 

xWjGo^^O + eGi^z)} =0. (2-8) 

G +eGi=0 

Because we wish to find the approximate symmetries to leading order, we need only 
solve the following leading-order determining equation: 

X^{G (z,z,zf^) + eG 1 (z,z)} =0(e). (2-9) 

Go=0 

By solving the approximate determining equation, (|2-9|) . we obtain the infinitesimal 
generator X admitted by the system (|2-3|) in the leading-order approximation. As 
shown below in some examples, the approximate solutions of the determining equa- 
tion Q2-9p give various non-trivial approximate symmetries even if the system ()2-3|) 
admits only a few exact, trivial symmetries. It is by virtue of this fact that the 
present method is useful in obtaining approximate solutions. Using the approximate 
infinitesimal generator X, we construct a group- invaritant solution of the system 
(EH). 

The group-invaritant solution, z = z(e, t) , which satisfies 

X{z-z(e,t)} =0, (2-10) 

Z — Z{£,t)=0 
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is obtained by solving the so-called Lie equation 



d £ z 



£{t,z,z)d t z + n z (t,z,z). 



(2-11) 



Because our symmetry (|2-5j) always has a non-vanishing component in the e direc- 
tion, the Lie equation becomes a system of differential equations with respect to 
the perturbation parameter e as described above. Noting this fact, we refer to the 
Lie equation 1)2-11)1 as the renormalization group equation. Solving the renormaliza- 
tion group equation by adopting solutions of the unperturbed system as boundary 
condition, i.e. 



we obtain an approximate solution of the system 1)2-3)1 for non-vanishing e, 



As the solution thus obtained is invariant with respect to the Lie group, we call it 
the invariant renormalized solution. 

Here, some detailed discussion is necessary to regorously construct our remor- 
malization group method with the Lie symmetry. In the above-described general 
procedure, we obtain the approximate symmetry only to leading order in e, that is, 
e . Nevertheless, the solution invariant with respect to the leading-order Lie sym- 
metry yields an approximate solution up to 0(e) of the system 1)2-3)1 . Is this actually 
a solution of the given differential equations to order e 1 ? The following proposition 
answers this question. 

Proposition. Consider a system of ordinary differential equations of the form 



Assume that (12-141) approximately admits a Lie group transformation with an in- 
finitesimal generator X, i.e. 






(2-13) 




(2-14) 




(2-15) 



where X^ is the first prolonged infinitesimal generator. Then, its invariant solution 
is a solution of 1)2-14)1 to order e 1 
Proof. Suppose that 1)2-14)1 admits 



X = d £ + z, z)d t + v z (t, z, z)d z + V z (t, z, z)d s . 



(2-16) 



Then its invariant solution is a solution of 



X{z - z(e,t)} 



= 

z=z(e,t) 



(2-17) 




-£(t, z, z)d t z(e, t) + rj z (t, z, z),v(s, t) 



S(t,u,v)dt (2-18) 
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Thus, the invariant solution can be expressed as 

z(e,t) = z ( -°\t)+ [ {-£(t,z,z)d t z(e,t)+r 1 z (t,z,z)}de (2-19) 
Jo 

= z (0)_ j E ^(t,z^,z^)z^' - v z (t,z^,z^))de + 0(e 2 ) (2-20) 

= z® - e fat, z (0 \z^)z^' - r, z (t, *(°), f <°>)) + 0(e 2 ), (2-21) 

where z^ (t) and z^ (t) are the solutions of the unperturbed system. By substituting 
1)2-21(1 into (|2T4[) and taking account of (|2-15() . we can show that the solutions satisfy 
the given differential equations to order e 1 □ 

It is noted that this proposition is easily proved for a solution of arbitarary order 
in e by expanding (|2-19|) to that order. Owing to this proposition, it is not necessary 
to know a naive perturbed solution in order to find the proper Lie symmetry. This 
contrasts with the situation in the RG method of Shirkov and Kovalev, 5 ) 6 ) in which 
the naive perturbed solution is needed. 

§3. Examples 

Let us consider two examples. 
(1) a Harmonic oscillator equation 

Here, we apply the renormalization group method with the Lie symmetry to the 
equation for a harmonic oscillator as an example of singular perturbaiton problems. 
We consider the equation 

u" + u = -eu. (3-1) 
Then, defining z := u + iu', we have 

z' = -iz - + z) (3-2) 

Next, let 

X(t, z, z) = d £ + £(t, z, z)d t + n z (t, z, z)8 z + rf{t, z, z)3- z (3-3) 

be the infinitesimal generator of an approximate Lie symmetry admitted by (|3-2|) . 
Then, we write its first prolonged infinitesimal generator as 

X« (t, z, z, z', z') = d £ + z, z)d t + r] z (t, z, z)d z + rf{t, z, z)d- z 

+r ? 2 « (t, z, z, z', z')d z , + r? 2{1) (t, z, z, z', z')d- z , , (3-4) 

where 

Tf z, z, z', z') = V z t + z'inl - i t ) + z'r,l - z'\ z - z'z%, (3-5) 
V zil \t, z, z, z', z') = 4 + z'(4 - 6) + z'4 - z' 2 i- z - z'z'i z . (3-6) 
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Here, the subscript a (a = t,u,v) of ^and n represents the operation of d a . Then, 
the approximate determining equation for ()3-2j) is 



\z' + iz + e-i(z + z)\ =0(e 

{ 2 J z'+iz=0 



(3-7) 



This equation leads to 



+irf + U(z + z)=0. (3-8) 

Z' + IZ=0 Z 



The system (j3-8j) is linear in the unknown rf and £ with an inhomogeneous term 
due to the perturbation. Because we are interested in the effect of the perturbation, 
we seek a particular solution of the linear inhomogeneous system (|3-8|) . Expanding 
£ and n u in powers of t, z and z as 

= £ H jk lt j z k z l , (3-9) 
rf(M^) = ^ Vjkit j z k z l . (3-10) 

j,k,l>0 

and substituting these expressions into (|3-8|) . we obtain a particular solution of lowest 
order, 

S(t,z,z) = 0, (3-11) 
77 z (t,z,z)=-ifc + iz-iz. (3-12) 

Then, the infinitesimal generator admitted approximately by (|3-2p is 

X(t, u, v) = d e + ( —-itz + \z — -z] d z + (\itz — \z + -z] cV (3-13) 



The invariant solution is derived from the corresponding Lie equation or renormal- 
ization group equation, and we have 

d e z = ~itz + - A z-^z. (3-14) 

This renormalization group equation is consistent with the naive expanded solution 
of (l3~21l . 

z(e,t) = Ae'^+V+e (-^iAte'^ + ^Ae'^ - ±A e <(*+')) + 0(e 2 ).(3-15) 
Then, defining z := Ae~^ t+e \ Ea. ^T^i reads 

z(e,t) = z + e (--itz Q + -z + -z j , (3-16) 
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which is the Euler scheme difference equation corresponding to the renormalization 
equation (|3-14j) . Using the unperturbed solution as the boundary condition, i.e. 

z(e = 0,t) = Ae- i{t+e \ (3-17) 

we obtain the invariant renormalized solution 

u(e,t) = 2 ( t 2 - ~ ) " etiisin ( - ( t 2 - - ) " e ] (2t cos(i + 9) - sin(i + 0)) 




+e4i?cos - r - - £\sm(t + e), (3-18) 



where R = \A\. Although the expression (|3-18|) is somewhat complicated, it is 
easy to show that this invariant solution asymptotically approaches the renormalized 
solution derived using the conventional RG method. Since we are interested in the 
long-time behavior of solution [i.e. t ~ in this case] we set (i 2 — |) ~ t 2 and 

e* k 1 and ignore sin(i + 9) in the first term. Then, the solution (|3-18|) reduces to 
the conventional renormalized solution, 

u(e,t) =Rsin ((l + e^)t + 9j. (3-19) 

Adding the general solution of the determining equation ()3-8|) to the particu- 
lar solution given in Eas. (|3-ll|) and (|3-12j) . we can construct a simpler invariant 
renormalized solution. For example, choose a particular solution of (|3-8|) as 

= -~t, (3-20) 

rj z (t,u,v)=^(z-z). (3-21) 

Then, the Lie equation becomes 

d E z = ~td t z + -(z-z). (3-22) 

Thus, we have the simpler invariant renormalized solution 

u(e,t) = Rsm(eU + 9), (3-23) 

which also is identical with the conventional renormalized solution to order e 1 . These 
two invariant renormalized solutions are depicted in Fig. 1 along with the exact 
solution, the naive perturbation solution and the conventional renormalized solution. 
We observe good agreement among these solutions even for e = 0.2, except for the 
naive perturbation solution. 

(2) Rayleigh equation 
We next consider the Rayleigh equation, 

u" + u = e ( u' - -u' 3 ] , (3-24) 



Asymptotic Solutions Invariant under the Renormalization Group 



7 




Fig. 1. Comparison of perturbative solutions of Eg. 13-11 . (1) represents the naive perturbation 
solution of order e, (2) the exact solution, (3) and (4) the invariant renormalizad solutions l|3-18fl 
and 113-2311 . respectively, and (5) conventional renormalized solution, where 9 = 0, R = l,and 
e = 0.2 



which can be rewitten 

z* = -iz + e ( -z z + —z 3 - -\z\ 2 z + -\z\ 2 z - —z 3 ) . (3-25) 

V 2 2 24 8 1 1 8 1 1 24 J y ' 

In this case, the determining equation for the components (£, rf ,rj z ) of the infinites- 
imal generator is 



2 1 1 1,1.., 1, ,o 1 , 

+in z + -z--z + —z 3 - -\z\ 2 z + -\zfz -z 3 = 0. 3-26 

z'+iz=o 1 2 2 24 8 1 1 8 1 1 24 y J 



Following the same procedure as in the previous example, we obtain the following 
particular solution of (|3-26|) : 

£(t,u,v) = 0, (3-27) 

n z (t, z, z)= -iz + -tz + —iz 3 - —iz 3 - —i\z\ 2 z - -t\z\ 2 z. (3-28) 
' v ' ' ; 4 2 48 96 16 1 1 8 1 1 y J 

From this, we find the infinitesimal generator approximately admitted by Eg. ()3- 24(1 
to be 

X(t,u,v) = d £ + ( -iz + -tz + —iz 3 — —iz 3 — —i\z\ 2 z - -t\z\ 2 z) d z 
v ; V 4 2 4 8 96 16 1 1 8 1 1 / 

1 1 1 q 1 <i 1 , ,9 l M ,\ n ( 

--iz + -tz H iz 3 iz 3 H i\z\ 2 z t\z\ 2 z d- z . (3-29) 

4 2 96 48 16 1 1 8 1 1 J y 1 

The corresponding Lie equation or renormailization group equation reads 

d £ z = -iz + -tz + —iz 3 - —iz 3 - —i\z\ 2 z - -t\z\ 2 z. (3-30) 
£ 4 2 48 96 16 1 1 8 1 1 y J 
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This renormalization group equation is also consistent with the naive expanded so- 
lution of dSZ2U), 2) 

z = Ae^+V 

+^( 1 -T)^ ( " ) + l( 1 -T)^ +e) 

+ ^ ie -3i(t + e) + ^ ie si(t + e)^ (3 . 31) 

Defining z := Ae~ i(t+e \ Eq.(j33U) reads 

f 1 f 1 o 1 , 1 , ,2- 1 , ,2 1 

2 = zo + ej + -jtzo + — + — iz Q - -^Azq\ z - -t\z \ z V , (3-32) 

which is the Euler scheme difference equation corresponding to the renormalization 
equation (|3-3U|) . The renormalization group equation (|3-3U|) is approximated for large 
t by 

d e z=±tz(l-j\z\A. (3-33) 
Under the transformation of coordinates (z,z) i— > (A, a), with 

A:=\z\, a:=|Log(~), (3-34) 



Ea. (|3-33j) becomes 



or 



d £ A = 4(i-4M/. 

d F a = 0. 



d T A = 
d T a = 0, 



(3-35) 



(3-36) 



where r := et. This renormalization group equation is equivalent to that obtained 
with the conventional renormalization group method. 1 ) 

Because the determining equation of the Lie symmetry is linear in each compo- 
nent of the infinitesimal generator, it is easy to calculate higher-order corrections to 
the leading-order Lie symmetry. The determining equation to order e 1 is 



z +iz-e{ 2 ~ 2 ~ , 24 «. s i~i 8 i-i - 24* 

-e ( -z - -z + ^-z 3 - -\z\ 2 z + -\z\ 2 z- | =0(e 2 ) (3-37) 
V2 2 24 8 1 1 8 1 1 24 J v 1 y 1 

The solution of (|3-37|) gives the following approximate symmetry to order e 1 for large 
t: 

V* = \tz(l-\\zA + e\tiz(l + \\z\ 2 -±\zA. (3-38) 
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Then, we obtain the renormalization group equation to order e 2 , 

| = ^(i- jN 2 )t + 4«(i + jN 2 -^N 4 )*, (3-39) 

or 

§4. Concluding remarks 

In the method presented here, the renormalization equation appears as the main 
ingredient in the Lie group theory. More specifically, we find that the Lie equation 
for the Lie point group in the expanded space that includes the perturbation pa- 
rameter is the renormalization group equation. In the conventional RG method, it 
is necessary to calculate naive perturbation solutions of a given nonlinear system, 
which is often tedious, due to the nonlinearity of the system. To avoid this step, the 
proto-RG approach was proposed in the perturbation analysis of a perturbed sys- 
tem. 7 ) The present method does not require any pertubational analyses to determine 
approximate solutions of a perturbed system but, instead, approximate symmetries. 
Because the determinig equation is always a linear system for the symmetry, it is 
easier to calculate an approximate symmetry than to obtain an approximate solu- 
tion of the nonlinear system. Thus, the present method always remains within the 
framework of Lie group theory and completely frees us from the need to calculate 
naive perturbation solutions. 

It should be noted that the present method can be applied not only to continuous 
systems but also to discrete systems, owing to the general framework of the Lie group. 
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